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Abstract. For a supersonic Euler flow past a straight wedge whose vertex 
angle is less than the extreme angle, there exists a shock-front emanating from 
the wedge vertex, and the shock-front is usually strong especially when the 
vertex angle of the wedge is large. In this paper, we establish the L 1 well- 
posedness for two-dimensional steady supersonic Euler flows past a Lipschitz 
wedge whose boundary slope function has small total variation, when the total 
variation of the incoming flow is sufficiently small. In this case, the Lipschitz 
wedge perturbs the flow and the waves reflect after interacting with the strong 
shock-front or the wedge boundary. We first obtain the existence of solutions 
in BV when the incoming flow has small total variation by the wave front 
tracking method and then study the L 1 stability of the solutions. In particular, 
we incorporate the nonlinear waves generated from the wedge boundary to 
develop a Lyapunov functional between two solutions, which is equivalent to 
the L 1 norm, and prove that the functional decreases in the flow direction. 
Then the L 1 stability is established, so is the uniqueness of the solutions by 
the wave front tracking method. Finally, we show the uniqueness of solutions 
in a broader class, i.e. the class of viscosity solutions. 



1. Introduction 

For the Cauchy problem of a strictly hyperbolic system of conservation laws: 

(1.1) U t + F{U) a = 0, UeR n , 

(1.2) U\t=o = U(x), 

whose each characteristic field is either linearly degenerate or genuinely nonlinear, 
the existence of weak solutions to (|1.1|) - (|1.2|1 with small total variation was first 
proved by Glimm by a probabilistic algorithm, the Glimm scheme; and a de- 
terministic version of the Glimm scheme was developed by Liu |28| . Alternative 
methods for constructing solutions of the Cauchy problem were first introduced in 
[1(51 118| . based on wave front tracking. For the scaler equation, F is approximated 
by piecewise linear functions F v in Dafcrmos so that the approximate solutions 
are piecewise constants and all the interactions are determined by solving the Rie- 
mann problem. The method was generalized to the 2x2 case in DiPerna ^H] m 
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which piecewise constant approximate solutions are constructed so that the wave 
interactions can be determined by only solving the Riemann problem. Bressan |2j 
developed the wave front tracking method for nxn systems by overcoming the diffi- 
culty that the procedure used in |18| may yield an infinite number of discontinuities 
in finite time when n > 2; and the wave front tracking method was further simpli- 
fied later in Baiti-Jenssen pQ. Also see Bressan [5], Dafermos ^7], Holden-Risebro 
|21|. and LeFloch |22] for further references. 

Within the class of initial data U 6 L 1 n BV(R; R n ) with suitably small total 
variation, it was established that problem (|1 — (fl.^ft is well-posed in L 1 (i?;i?™) 
for the solutions generated by the wave front tracking algorithm. In particular, in 
Bressan- Colombo Bressan-Crasta-Piccoli [S], and Bressan-Liu-Yang 0, it was 
proved that the entropy solutions of 1)1.1(1 - 111.2(1 constitute a semigroup which is 
Lipschitz continuous with respect to time and initial data. Lewicka-Trivisa |26j 
obtained the L 1 well-posedness of solutions generated by the wave front tracking 
method for the Cauchy problem 1(1.1(1 1(1.2(1 with the initial data U being a small 
perturbation of a fixed Riemann problem ([/_,[/+) containing two large shocks, 
under the necessary stability condition (cf. [3 ESI EH; also see |23 ESI)- The L 1 
well-posedness in the class of viscosity solutions for the Cauchy problem has been 
also established (cf. and the references therein). 

In this paper, we are concerned with the L 1 well-posedness of a physical nonlin- 
ear problem of initial-boundary value type, which governs two-dimensional steady 
supersonic Euler flows past a curved wedge. More specifically, the two-dimensional 
steady supersonic Euler flows are generally governed by 



(1.3) 



(PU) X + {pv)y = 0, 

(PU 2 +P) X + (pUv)y = 0, 

(PUV ) X + (pV 2 + p)y = 0, 

{pu{E + p/p)) x + (pv(E + p/p)) y = 0, 



where (it, v) is the velocity, p the density, p the scalar pressure, and E = \{ u + 
v 2 ) + e(p,p) the total energy with e the internal energy (a given function of (p,p) 
defined through thermodynamical relationships). The other two thermodynamic 
variables are the temperature 8 and the entropy S. For an ideal gas, 

R 

(1.4) p = Rp8, e = c v 9, j = l-\ > 1, 



Cv 



and 



(1.5) P = p(ftS) = K pV/ c », e = -Z-pi-W* = 

7 — 1 7 — 1 

where R, k, and c v are all positive constants. 

If the flow is isentropic, i.e. S — const., then the pressure p is a function of the 
density p, p = p(p), and the flow is governed by the following simpler isentropic 
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Euler equations: 

f {PU) X + (pv)y = 0, 

(1.6) < (pu 2 +p) x + (puv)y = 0, 

[ (PUV) X + (pV 2 + Pjy = 0. 

For polytropic isentropic gases, by scaling, the pressure-density relationship can be 
expressed as 

(1.7) p(p)=pVl, 7>1- 
For the isothermal flow, 7=1. The quantity 

c = \Jp P (Pi S ) 

is defined as the sonic speed and, for polytropic gases, c = \pjpfp. 

System (|1.3H or i|1.6fl governing a supersonic flow (i.e., u 2 + v 2 > c 2 ) has all real 
eigenvalues and is hyperbolic, while system 1)1. 3|) or (|1.6J) governing a subsonic flow 
(i.e., u 2 + v 2 < (?) has complex eigenvalues and is elliptic- hyperbolic mixed and 
composite. 

The study of two-dimensional steady supersonic flows past a wedge can date back 
to the 1940s (cf. Courant-Friedrichs 53)- Local solutions around the wedge vertex 
were first constructed in Gu [201; Li [22], Schaeffer (2H|, and the references cited 
therein. Global potential solutions were constructed in various different setups in 
[1 H 1121 [TS1 131) when the wedge vertex angle is less than the critical angle. 

For the full Euler equations, when a wedge is straight and the wedge vertex 
angle is less than the critical angle, there exists a supersonic shock-front emanating 
from the wedge vertex so that the constant states on both sides of the shock- 
front are supersonic; the critical angle condition is necessary and sufficient for 
the existence of the supersonic shock (cf. Courant-Friedrichs ^E])- When the 
incoming flow is uniform, Chen-Zhang-Zhu ^U] first established the existence of 
global supersonic Euler flows, especially the nonlinear structural stability of the 
strong shock-front emanating from the wedge vertex under the BV perturbation 
of the wedge boundary. In this paper, we first show the existence of solutions to 
the above problem when the incoming flow is a BV perturbation of the uniform 
flow by the wave front tracking method, and then we establish the ^-stability of 
entropy solutions generated by this method. Based on these, we establish estimates 
on the uniformly Lipschitz semigroup V generated by the wave front tracking limit 
and prove the uniqueness of solutions by means of local integral estimates within a 
broader class of solutions, i.e. the class of viscosity solutions. 

One of the main new ingredients in our approach here is to develop techniques to 
handle the boundary difficulty, in comparison with the earlier works on the Cauchy 
problem. For the L 1 stability of solutions of the Cauchy problem, the decrease 
of the Lyapunov functional was achieved by essentially using the cancellation of 
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the distances on both sides of waves. However, for our wedge problem that is 
the problem of initial-boundary value type, there is no such cancellation near the 
boundary, since only one-side is possible near the wedge boundary. In order to 
overcome this difficulty, we employ the exact feature of the boundary condition to 
obtain an estimate and refine the functional based on this estimate. In particular, 
since the flow of two solutions near the boundary must be parallel, we identify the 
relation between these two states, which is desirable for redesigning the functional 
to ensure the decreasing of the functional in the flow direction. 

For concreteness, as in Chen-Zhang-Zhu JTQJ , we will analyze the problem in the 
region below the lower side T of the wedge for the Euler flows for U — (u,v,p,p) 
governed by system (II .3|) and U = (u,v,p) by the case above the wedge can 

be handled in the same fashion. Then we have 

(i) There exists a Lipschitz function g G Lip(R + ) with g' G BV(R + ), g'(0+) = 
0, and g(0) = such that 

^■■={{x,y) : y < g(x), x > 0}, V := {(x,y) : y = g(x), x > 0}, 

and n(x±) = fcl ( x± ^^ is the outer normal vector to T at the point x± 
(see Fig. EJ; 




O 

Figure 1 . Supersonic flow past a curved wedge 

(ii) The upstream flow U = U(y) = (u, v,p, p)(y) at x = satisfies 
u(y)>0, v(y)>0, u(y) 2 + v(y) 2 > c(y) 2 := 7p(y)/ p(y), 
and 

< arctan(s(y)/«(y)) < u) cr it, 

where uj C rit is the critical vertex angle so that there is a supersonic shock- 
front emanating from the wedge vertex. 

With this setup, the wedge problem can be formulated into the following problem 
of initial-boundary value type for system (|f .3(1 or (|I.6(I : 
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Cauchy Condition: 

(1.8) U\ x = = U(y); 
Boundary Condition: 

(1.9) (u,v) ■ n = on T. 

This paper is organized as follows. In Section 2, we discuss the basic proper- 
ties for the adiabatic Euler equations and related nonlinear waves. In Section 3, 
we discuss the wave front tracking method and define the interaction potential Q, 
and then we prove the existence of entropy solutions to the initial-boundary value 
problem. In Section 4, we construct the Lyapunov functional $ between two solu- 
tions to incorporate the nonlinear waves generated by the wedge boundary vertices, 
which is equivalent to the L 1 distance between these two solutions. In Section 5, 
we prove the decrease of $ in the flow direction, which implies the L 1 stability 
of the solutions. In Section 6, we prove the existence of the semigroup generated 
by the wave front tracking method and establish some estimates on the uniformly 
Lipschitz semigroup S generated by the wave front tracking limit. In Section 7, 
we prove the uniqueness of solutions by means of local integral estimates within a 
broader class of solutions. 



2. Euler Equations and Nonlinear Waves 

In this section, we review some basic properties of the adiabatic Euler equations 
(|1.3|) and related nonlinear waves, which will be used in the subsequent sections. 
The Euler equations for steady supersonic flows can be written in the following 
conservation form: 

(2.1) W(U) X + H(U) y = 0, U ={u,v,p,p) T , 

with 



W(U) = (pu,pu 2 +p,puv,pu(h-\ j) T ,H(U) = (pv,puv,pv 2 +p,pv(h+ )) T , 



and h — , 2^)^ ■ The eigenvalues of system (|2.1|) are 



2.2 I, = ^ — -^r— 2 , 3 = 1, 4; U = v u, i = 2, 3, 

where c 2 — jp/p. If the flow is supersonic (i.e. u 2 + v 2 > c 2 ), system (|2.1|) is 
hyperbolic; and, in particular, when u > c, system l|2.1f) has the four corresponding 
linearly independent eigenvectors: 

(2.3) v 3 ■ = K 3 (-l J ,l,p{l J u-v) lP {l 3 u-v)/c 2 ) T , i = l,4, 
r 2 - (u,v,0,0) T , r 3 = (0,0,0, p) T , 
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where chosen so that Tj ■ Vlj — 1 since the jth-characteristic fields are 

genuinely nonlinear, j = 1,4. Note that the second and third characteristic fields 
are always linearly degenerate: Tj ■ VAj = 0, j = 2, 3. 

Definition 2.1 (Entropy Solutions). A BV function U = U(x,y) is called an 
entropy solution of the initial-boundary value problem <|2.1[1 and <|1.8[l - <|1.9|) provided 
that 

(i) U is a weak solution of l)2.1|l and satisfies 

U\ x =o = U(y), (u,v) ■ n\ y=g ^ =0 in the trace sense; 

(ii) U satisfies the entropy inequality, i.e. the steady Clausius inequality: 

(2.4) (puS) x + (pvS) y > 

in the sense of distributions in fi including the wedge boundary. 

We now discuss the wave curves in the phase space. The contact Hugoniot curves 
Ci(Uo) through Uq are 

(2.5) Ci(U ) : p = p , w = v/u = v /u , « = 2,3, 

which describe compressible vortex sheets. We remark that, although the two 
contact discontinuities coincide as a single vortex sheet in the physical xy-plane, 
it requires two independent parameters to describe them in the phase space U = 
(u,v,p,p) since there are two linearly independent eigenvectors corresponding to 
the repeated eigenvalues A2 = A3 = v/u of the two linearly degenerate fields. 

Moreover, the rarefaction wave curves Rj(Uo) in the phase space through Uq are 



R (Uq) '■ dp = c 2 dp, du ~ —Xjdv, p(XjU — v)dv = dp for p < po, j = 1, 4. 



(2.6) 

The Rankine-Hugoniot conditions for (|2.1|l are 
(2.7) s[W(u)] = [H(u)}, 

where s is the propagation speed of the discontinuity. Then 



u vo + (-iycoy/i4 + v% -cl . . 
s = Sj := 2 — =2 '3 = 1)4, a = Oi = vo/u ,i = 2,3, 

Uq — Cq 

where c 2 Q = $~B- and b = ^ - 

u b po 1 2 2 po 

Plugging Si,i = 2,3, into l|2.7() . we obtain the z^-contact Hugoniot curves 
Ci(Uo),i — 2,3, as defined in (|2.5|l : while plugging Sj,j = 1,4, into l|2.7|l . we 
obtain the j t,l -Hugoniot curve Sj(Uo) through Uq: 

c 2 

Sj(U ) : [p] = T-[p], M = Sj[v], p (sjU -v )[v] = [p] for p > p , j = 1, 4. 
oo 
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The half curves of Sj(Uo) for p > po, denoted by Sj~(Uo),j = 1,4, in the phase 
space are called the shock curves on which any state with Uq forms a shock in the 
x — y plane satisfying the entropy condition as explain in Lemma 2.1 below. 

Note that Sj(Uo) contacts with R~ (Uo) at Uo up to second-order. 

As indicated in we have 

Lemma 2.1. If U is a piecewise smooth solution, then, on the shock wave, the 
entropy inequality H2.4(l in Definition 2.1 is equivalent to any of the following: 

(i) The physical entropy condition: the density increases across the shock in the 
flow direction, 

(2.8) P front < Pback', 

(ii) The Lax entropy condition: on the j th -shock with the shock speed o~j, 

(2.9) \j(back) < Sj < X^front), j = l,4, 

(2.10) si < X 2 , 3 (back), X2^(front) < s 4 . 

The following properties and related estimates of wave interactions in Lemmas 
2.2-2.8 have been obtained in Chen-Zhang-Zhu ^Oj- We list them below for subse- 
quent use in this paper. 

2.1. Riemann problems and Riemann solutions. We start with Riemann 
problems and their solutions. 

Lateral Riemann problem. The simplest case of problem (|2.1f) and p.8|l - p.9|l 
is g = 0. It can be shown that, if g = 0, then problem l|2.1|l admits an entropy 
solution that consists of a constant state U- and a constant state [/+, with U+ = 
(u+, 0,p+, p+) and u + > c + > in the subdomain of separated by a straight 
shock emanating from the vertex. That is to say that the state ahead of the shock- 
front is U-, while the state behind the shock- front is U+. When the angle between 
the flow direction of the front state and the wedge boundary at a boundary vertex 
is larger than it, then an entropy solution contains a rarefaction wave that separates 
the front state from the back state. 

Riemann problem involving only weak waves. Consider the following initial value 
problem: 

W(U) x + H(U) v = 0, 
U a , 

U b , y < yo, 
where Ub and U a are constant states. 

Lemma 2.2. There exists e > such that, for any states U a ,Ub £ O e (U+) or 
U a ,Ub € O e (U-), problem (|2.11|l admits a unique admissible solution consisting of 
four elementary waves. 



(2.11) { p., = u _ / U a , y > 2/0, 
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Riemann problem involving a strong 1-shock. For simplicity, we use notation 
{Ub, U a } = (ai, «2) &3, O-a) to denote the solution to the Riemann problem, where 
cti is the strength of the ith elementary wave. For any U G 5i ({/_), we also use 
{U-, U} = (a, 0, 0, 0) to denote the 1-shock that connects U- and U with speed a. 
Then we have 

Lemma 2.3. Let {[/_, U+} = (oo, 0, 0, 0), p + > /9_, and 7 > 1. Then 
o~ < 0, u + < u_ < (1 + 1/7)14+, 

and 

det(V £/ F(?7 + ) - a VuW(U + )) > 0. 

Furthermore, there exists a neighborhood O e {U+) ofU+ and a neighborhood O e (U-) 
of U- such that Uq G O s (UJ) and the shock polar Si(Uq) D O e (£/+) can be param- 
eterized by the shock speed a as a — > G{Uq,o) with G G C 2 near (f7_,cro) and 
G(U-,<t ) = U+. 

2.2. Estimates on wave interactions and reflections. We have 

Lemma 2.4 (Estimates on weak wave interactions). Suppose that Ub,U m ,U a G 
O e (U+), or Ub, U m , U a G O s (U-), are three states with {Ub, U rn } = (ai, a<x, CK3, a^), 
{U m ,U a } = (/3i,/3 2 ,/33, AO, and {U b , U a } = (71,72,73,74)- Then 

7 i = a i +/3 i + 0(l)A(a > /3), 

w/iere A(a,/3) = (|a 4 | + |a 3 | + |a 2 |)|/3i| + M(|/? 2 | + |As|) + Ej=i, 4 A j( a > Z 3 ) 

■ ,v ; [ otherwise. 

Denote {Ck{a k , b k )}kLo the points {(afc, £>jt)}j£L in the xy-plane with dfc+i > 
a k > 0. Set 

/ &fc+i — &fc \ 

Wfe^+i = arctan , u k = ujk,k+i - ^k-i,k, U-1,0 = 0, 



(2.12) O fc+ i = {(x,y) : x <E [a k ,a k +i), y < b k + (x - a k ) taxi(u k ,k+i)} , 
T/c+i = {(x, y) : x G [a fc , Ofc+i), y = fofe + (x - a fc ) tan(wft,fc+i)}, 

and the outer normal vector to T k '- 

/ 010 % (6fc — i»fc+i,Ofc+i — Ofc) t ■ 1 \ { w 

(2.13) n fe+ i = — = (-sm(w fc , fc+1 ),cos(a; i!>fe+ i)). 

V(bfc+i - 6 fe ) J + (ofc+i - afe)^ 
Then we consider the initial-boundary value problem with U_ a constant state: 

GU in0 fc +i, 
U\ x=ak =U, 

(u,v) ■ nfe+i =0 on T k+ i. 
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Lemma 2.5 (Estimate on the boundary perturbation of the strong shock). For 
e > sufficiently small, there exists e — s(e) < e so that G(Oi(o~o)) C O e (U+) 
and, when \0 k \ < s, the equation G(a) ■ (nfc,0, 0) = admits a unique solution 
o~k G Oi(aa). Moreover, we have 

(2.14) (Tfc+l = cr fe + K bs k , 

where \K bs \ is bounded. 

Lemma 2.6 (Estimate on the boundary perturbation of weak waves). Let = 
(life, Vk,Pk, pk) be the state near the boundary with (u^, Vk)-tik = 0. Then there exists 
Uk+i such that Uk+i} = (Si, 0, 0, 0) and (uj^+i, ffe+i) • n^+i = 0. Furthermore, 

Si = K bQ k , 

where K b $ is bounded. 

Lemma 2.7 (Estimates on the reflection of weak waves on the boundary). Let 

{U b ,U k } = (0,02,03,0:4) and (uk,Vk) ■ n/. = 0. Then there exists Uk+i such that 
{U b ,U k +i} = (Si, 0,0,0) and (u k+1 ,v k+ i) ■ n k = 0. Furthermore, 

Si = K b4 a 4 + K b3 a 3 + K b2 a 2 , 

where K b4 , K b 3, K b2 , and K b o are C 2 — functions of (04, 03, a 2 ,0i, 0k', U b ) satisfying 

Kb4\{0 k =ci A =ot3=oi2=Pi=Q,U h =U + } = 1j 

K b2 \[ 0k=a4=a3=a2= p 1=Oi u b =u + } = K b3 \[ 0k=ai=a3=a2= p 1=o ^ Ub=u+ j = 0. 

Lemma 2.8 (Estimates on the interaction between the strong shock and weak 
waves from above). Let U m ,U a G O e (U+) with {G(U b ,cr),U m } = (0,0,0,0) and 
{U m ,U a } = (01,0, 0, 0). Then there exists a unique (a', S 2 , S3, 84) such that the 
Riemann problem (|2.11(l with U b G O e (U~) admits an admissible solution consisting 
of a strong I -shock, two contact discontinuities of strengths S 2 and S3, and a weak 
A-wave of strength S 4 : 

{U b ,U a } = (a', 8 2 ,S 3 ,S 4 ). 

Moreover, 

a' = a + K s iPi, 8 2 = K s2 (3i, 83 = K s3 (3i, 84 = K s4 (3i, 

where \K s4 \ < 1, and \K s i\ + \K s2 \ + \K s3 \ is bounded. Furthermore, we have 

, 01 _, ~ Z 4 + - go , ctqM-Q - u + l 4+ P 

(2-15) K s4 \- = — — - <1. 

ti+ — a o- U-Q + u + l 4+ F 

Lemma 2.9 (Estimates on the interaction between the strong shock and weak 
waves from below). Let U m , U b G O e (U-) and U a G O e (U+) with 

{U b ,U m } = (01,02,03,04), {U m ,U a } = (a, f3 2 , fa, (3 4 ). 
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Then there exists a unique (a' ,82,63,64) such that the Riemann problem (|2.11(l 
admits an admissible solution consisting of a strong 1-shock, two contact disconti- 
nuities of strengths 82 and 83, and a weak 4-wave of strength 84: 

{U b ,U a } = (a', 82,83,64). 

Moreover, 

4 4 
u' = a + Kuoii + 0(1) A, 82 = p 2 + K * a i + °W A ' 

1=1 i=l 
4 4 

6 3 = I33 + Y, K ^ + °W A ' S * = ^ + 12 Ka ^ + W A ' 

i=l i=l 

where \Kji\,i,j — 1,..,4, are bounded and A = X)j=i 2 3 4 j—% 3 4 \ a i@j\- Further- 
more, we can write the estimates in a more precise fashion: 

4444 
0-' = a + 22 Kudu 8 2 = fh + ^ K 2i a " 63 = f3 3 +}^ K 3i ai, 84 = fa + 2J K 4l a l , 

i—1 i—1 i—1 i—1 

where YhJ=i < M f or some M > °- 

Proof. We first consider the interaction between (Ub,U m ) = (a±, 0,2, 0.3, C14) and 
(U m , G(U m , a)) — (a, 0, 0, 0) to find that the solution is the perturbation of the un- 
perturbed states of the strong shock. From Lemma lOl we know that U = G(Uo,cr) 
near (U-,o~o) with G G C 2 , which implies that U = (u,v,p,p) depends contin- 
uously on the state Uq = (uq, vo,po, po)- The perturbation near U— is equiva- 
lent to 53j=i 0(l)ai. Then the interaction estimate between (Ub,G(U m ,a)) and 
(G(U m , a), U a ) follows from Lemma |2~31 □ 

3. Wave front tracking method and existence of entropy solutions 

The basic idea of the wave front tracking method is to construct approximate so- 
lutions within a class of piecewise constant functions: First, approximate the initial 
data by a piecewise constant function and solve the resulting Riemann problems 
and replace the rarefaction waves by step functions with many small discontinu- 
ities; then track the outgoing fronts until the first time when two waves interact 
which are determined by a new Riemann problem; and finally design a simplified 
Riemann solver so that the number of wave fronts is finite for all x > in the flow 
direction. 

3.1. The Riemann solvers. As mentioned in Section 2, the solution to the Rie- 
mann problem (Ub,U a ) is a self-similar solution given by at most five states sep- 
arated by shocks or rarefaction waves. More precisely, there exists C 2 curves 
a — > iJj(oi)(u) parameterized by arc length such that 

U b = ^4(0^4) o . . . o ipxta^iUa) 
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for some a = (ai, . . . , 0C4) and Ui = ipi(ai) o . . .0 ipi(ai)(U a )- When on is positive 
(negative), states Ui-i and Ui are separated by an i-rarefaction (i-shock) wave, so 
we call at the strength of the z-wave. 

For given initial data U, let Lf £ , e > 0, be a sequence of piecewise constant 
functions approximating U in the L 1 norm, and the wedge boundary is also ap- 
proximated as in Section 2. Let N e be the total number of discontinuities in the 
function XJ e and the tangential angle functions of the wedge boundary. Choose a 
parameter S £ > controlling the maximum strength of rarefaction fronts, and A 
(strictly larger than all the characteristic speeds of H2.lt ) that is the speed of non- 
physical waves generated whenever the simplified method is used. The strength of 
the non-physical waves is the error due to the simplified Riemann solver. 

A. Accurate Riemann solver: The accurate Riemann solver is just the so- 
lution to the Riemann problem (as in Section 2), except every rarefaction wave 
(w, Ri(w)(a)) is approximated by a piecewise constant rarefaction fan. 

B. Simplified Riemann solver: For the weak waves, it is exactly the same as 
in pQ. When a weak wave interacts with the large shock, the simplified Riemann 
solver is that we ignore the strength of the weak wave, keep the strength of the 
strong shock, and put the error in the non-physical wave as follows: 

Case 1 (A weak wave (U-,Ui) hits the large shock (Ui,U+) from below): We 
solve the Riemann problem (U-, U+) in the following way: 

U- for y/x < A(U U U+), 

Ui for A(f/i, U+) < y/x < A, 

U+ for y/x > A, 

where A(t/i, U+) is the speed of the strong shock; 

Case 2 (A weak wave (U2,U+) hits the large shock (U-,1/2) from above): We 
solve the Riemann problem (f7_, £/+) in the following way: 

V- for y/x < A(U-,U 2 ), 

U 2 for A(U-,U 2 ) < y/x < A, 

U+ for y/x > A, 

where A(U-, U2) is the speed of the strong shock. 

3.2. The algorithm to construct the approximate solutions. Given e, we 
construct the approximate solution U £ (x,y) as follows. When x = 0, all the Rie- 
mann problems in JJ e are solved by the accurate Riemann solver. By slightly 
perturbing the speed of a wave, we can assume that, at any time, we have at most 
one collision involving only two incoming fronts. Let fj, e be a fixed small parameter 
with fi £ — + 0, as e — > 0, which will be specified later. For simplicity of notation, we 
will drop the index i in a* and do not distinguish between «j and a when there is 
no ambiguity from now on; the same for f3; also we use the same notation a as a 
wave and its strength as before. 



12 GUI-QIANG CHEN TIAN-HONG LI 

Case 1 (There is a collision between two weak waves with strengths a and (3 at 
some x > 0, respectively): The Riemann problem generated by this interaction is 
solved as follows: 

• If | a/9 1 > /ti e and the two waves are physical, then we use the accurate 
solver; 

• If \a(3\ < /u e and the two waves are physical, or one wave is non-physical, 
then we use the simplified Riemann solver. 

Case 2 (There is a collision between the large shock and one weak wave a at some 
x > 0): The Reimann problem generated by this interaction is solved as follows: 

• If \a\ > He and the weak wave is physical, then we use the accurate solver; 

• If |a| < [i e and the weak wave is physical, or this wave is non-physical, then 
we use the simplified Riemann solver. 

Case 3 (The wave hits the boundary or the boundary perturbs the flow): We use 
the accurate Riemann solver to solve the lateral Riemann problem. 

3.3. Glimm's functional and interaction potential. We now develop the Glimm- 
type functional and interaction potential for the initial-boundary value problem 
by carefully incorporating additional nonlinear waves generated from the wedge 
boundary vortices. 

Definition 3.1 (Approaching waves), (i) We say that two weak fronts a and (3, 
located at points x a < xp and belonging to the characteristic families i a ,ip € 
{1, . . . ,4} respectively, approach each other if the following two conditions hold: 

• x a and X/3 both lay in one of the two intervals into which R is partitioned 
by the location of the large 1-shock, i.e. the waves both belong to f2_ or£l + ; 

• Either i a = ip and one of them is a shock, or i a > ip. 

In this case we write (a, (3) G A. 

(ii) We say that a weak wave a of the characteristic family i a is approaching the 
large 1-shock if either a G f2_ and i a G {1,2,3,4}, or a G fl+ and i a — 1. We 
then write a G A\. 

(iii) We say that a weak wave a of the characteristic family i a is approaching 
the boundary if a £ f2_|_ and i a — 4. We then write a G Ab- 

For a weak wave a of i-family, we define its weighted strength as 

(a if a G 
a ~ \ fc_a if a G 

where k- = 2 max l-Kiij for coefficients Ku in Lemma \l. 91 

l<i<4,2<i<4 
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Definition 3.2. The wave interaction potential Q(x) is 

(3.1) Q(x) = c* ]T Kh\ + K * E E E l fcl 

(a,(i)eA aeAt 0£A b a k >x 

(3.2) = Q A + Qi + Qb + Q w , 
where K* e (K s4 , 1) and K b0 > K b $. 

Definition 3.3. The total (weighted) strength of weak waves in U e (x, •) is defined 
by 

v(x) = j2k\. 

a 

The Glimm-type functional is defined by 

(3.3) T{x) = V(x) + kQ(x) + \U*(x) -U+\ + \U.(x) -Uq\, 

where n > is a constant to be specified later, the vectors U*(x) and U*(x) are the 
above and below states of the large shock respectively at "time" x, and Uq and Uq 
are the right and left states of the large shock at x = 0, respectively. 

Note that V, Q, and T are constant between any pair of subsequent interaction 
times. On the other hand, we can show that, across an interaction "time" x, both 
Q and T decrease. 

Proposition 3.1. IfTV(U(-))+TV(g'(-)) is sufficiently small, then, for any x > 0, 
V(x) is sufficiently small and TV(U £ (x, •)) is uniformly bounded in e > 0. 

Proof. Define 

AT{x) = T{x+) -T{x~), 

where x + and x~ are the "times" right after and right before the interaction "time" , 
respectively. 

Case 1 (Weak waves a and [3 interact): Then U*{x) and U*(x) do not change 
across this interaction time. Thus, 

T(x + )-T(x-) = V(x+)-V{x-) + k{Q{x+)-Q{x-)) 

< MilU^I + K{-C*\b a bp\ + C*\b a b fi \V{x-) + M \b a b fi \), 

for some constants Mo and Mi independent of e. 

Case 2 (Weak wave a of 4-family interacts with the boundary): 

AJ r (x) = K b4 a -a + k (C*K b4 V(x~) a + K*K b4 a - a) . 

Case 3 (New 1-wave a produced by the boundary): 

AF{x) = K w k + k (c*K bo k V(x-) + K*K bo k - Kbo0k) ■ 

The next two cases when U*(x) and U*(x) change across this interaction "time". 
Then 
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Case 4 (Weak wave a of i-family interacts with the strong shock from below): 
AJ r (x) = V(x+)-V(x-) + \U*(x+) -U*(x~)\ 

+ \U*(x+) - U*(x+)\ + K(Q(x+) - Q(x~)) 
J2 K ri e a -b a + K(C* J2 K Jl V(x-)a-b a + K il a); 

3=1,2,3,4 3=2,3,4 

Case 5 (Weak wave a of 1-family interacts with the strong shock from above): 
AT{x) = K si a-a + K(C* K S3 V( X -) a - K*a + K s4 a) . 

.3 = 1,2,3,4 J=2,3,4 

In these cases, K S 4 < K* < 1, b a > 2max.{Kji}\a\ due to the choice of the weight 
k-, and C* > M > is a constant that is not small. 
We now prove 

V(x) < 1 for all x > 0. 
Case 1 (x > is the first interaction "time" x\): Since V(x^) = V(0) < 
TV{U{-)) < 1 and J27=o^k < TV{g'{-)) < 1 for all Cases 1-5, we find that, 
when k is larger enough and (j, s is sufficiently small, 

Af(xi) < 0, i.e. F( x t) < H x i) = H )- 

Therefore, 

V(x+) < T{x+) < J-(O) < V(0) + «Q(0) 

OO 

= y(o) + K (c*^ 2 (o) + T/(o) + i^^0 fc ) 

x=0 

OO 

< C(V(O) + J20k) «1. 

Case 2 (V(cc^) <C 1 and J r {x'^ ) < J r {x'^) for any fc < n): Then, for the next 
interaction "time" x„, similarly to Case 1, we also have 

AT(x n ) < 0, i.e. T(x+) < T{x-) = T(x+_ 1 ). 

Thus, we have 

V(x+) + \U*(x+) - U+\ + \U*{x+) - U Q -\ 

< no < n*n) = <■■■< ho) = v(o) + kq(o) 

oc 

= V(0) + K (C*V 2 (Q) + V(0) +KmY, x) 



x=0 



<C(V(O) + 5^0 fc ) «1. 



x=0 



Therefore, V(x) <C 1 is proved for all x, since C is independent of x. Then 
(3.4) TV{U(x, •)} « V(a;) + |t/*(z) - + |(7„(x) - f/ "| + K| = 0(1). 

□ 
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Lemma 3.2. For any sufficiently small fixed e > 0, the number of wave fronts in 
U E (x,y) is finite so that the approximate solutions U 6 (x,y) are defined for all x. 

Proof. Recall that the total interaction potential Q(x) is constant except decreasing 
when it crosses an interaction time. From Cases 1-5 in Proposition 13.11 we have 
known that V(t) <C 1. Therefore, we can find some c S (0, 1) so that 
(3.5) 

AQ(x) = Q(x+) - Q{x~) 

!— c\b a hp\ if both waves a and /3 are weak, 

— c|6 a | if weak wave a hits the strong shock or the boundary, 

— c|0fe | if the angle of the boundary changes. 

The following argument is similar to that in pQ: Q decreases for each case and Q(Q) 

is bounded; When the interaction potential between the incoming waves is greater 

than jU e , Q decreases by at least c/i e in these interactions, by the bound in l|3.5[) : 

Following the wave front tracking, new physical waves can be only generated by 

this kind of interactions, which implies that the number of the waves is finite; Since 

non-physical waves are produced only when physical waves interact, the number 

of non-physical waves is also finite; and, since two waves can interact only once, 

the number of interactions is also finite. Therefore, the approximate solutions are 

defined for all x > 0. □ 

Similar to £Q, we have following lemma. 

Lemma 3.3. The total strength of all non-physical waves at any x is of the order 
0(l)(5 c +n e ). 

Following the framework of the wave front tracking in 0^ and Lemmas I3.1H3.2I 
we obtain the existence of global entropy solutions to (|1.3(l and (|1.8(l - (|1.9(l . 

Theorem 3.1. IfTV(U(-)) + TV(g'(-)) is sufficiently small, then there exists a 
global entropy solution in BV of problem (|1.3|l and p. 8(1 - 1|1.9JI of initial-boundary 
value type in the sense of Definition 2.1. 

4. The Lyapunov functional 

We now follow the approach of EH EO] to construct the Lyapunov functional 
&(U,V) by incorporating additional new waves generated from the wedge boundary 
vortices, which is equivalent to the L^-distance: 

C^\\U(x, ■) - V(x, -)|| L i < *{U, V) < C 2 \\U(x, •) - V(x, -)\\ L i, 
^(U(x 2 r),V{x 2l -))~^{U{x 1 ,-),V{x 1 ,-)) < C 3 e(x 2 -x 1 ) for any x 2 > x x > 0, 
for some constants Ci, i = 1, 2, 3, where U and V are two approximate solutions ob- 
tained by the wave front tracking, and the small parameter e controls the following 
three types of errors: 
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• Errors in the approximation of initial data and boundary; 

• Errors in the speeds of shock and rarefaction fronts; 

• The maximum strength of rarefaction fronts; 

• The total strength of all non-physical waves. 

When x is fixed, for each y, the connection U(y) with V(y) always moves along 
Hugoniot curves Si, C%, C3, and S4 in the phase space. We call Pi(y) the strength 
of the z-th discontinuity wave, which is determined by U(y) and V(y) as follows: 

• If U(y) and V(y) are both in £7_ or in then start from U(y) moving 
along Hugoniot curves and end at V(y); 

• If U(y) is in V(y) is in Q + , then also start from U(y) moving along 
Hugoniot curves and end at V{y). 

• If V(y) is in J7_ and U(y) is in then start from V(y) moving along 
Hugoniot curves and end at U(y). 

Now we define the weighted L 1 strength: 

c iPi(y) if U(y) and V(y) are both in f2_, 

(4.1) = { cTPiiv) if U(y) and V{y) are in different domains, 

c iPi(y) if U(y) and V(y) are both in fi + , 



where the constants q, c™, and cf are to be determined later. Then we define the 
Lyapunov functional: 



4 .g( x ) 

(4-2) $(U,V)=J2 \Qi(v)\Wi(y)dy 

A — 1 J —OC 



with 



(4.3) Wi{y) = 1 + KiAi(y) + k 2 (Q(U) + Q(V)). 

Here Ki and K2 are two constants to be defined later, Q is the total wave interaction 
potential defined in Definition 13. 21 Ai(y) is the total strength of waves in U and V 
which approach the i-wave qi(y) defined by 



(4.4) Ai(y) = Bi(y) + Di(y) + 



Ci(y) if qi (y) is small, 

Fi(y) if i = 1 and <7i(y) = B is large, 
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where the "small" or "large" describes the waves that connect the states in the 
same or in the distinct domains 0~ and f2 + , respectively, and 



Bi(y) 



Ci(y) = 



E + E 

aej(U)Uj-(V) aej(U)Uj(V) 
Vol <V ,i<k a <4 y a >y,l<fca <i 



J2 ^(.u) +J2 c S j(v) )\a\ i£qi(y)<0, 



V ya<y,k 0i =i y a >y ,k a =i / 

fm = ( E + E 



aej(u)uj-(v) aej(u)uj(v) 

yo,<y,fc a — 1 Va^>y ik a — 1 

both states joined by a both states joined by a 
are located in 0_ are located in f2_|_ 



a . 



For fixed x, J = J(U) U J{V) is the set of all weak waves in U and V, a is the 
strength of wave a € J , located at point y a and belonging to the characteristic 
family k a . 



Di(y) = 


U, V are both in £!_ 


U, V are in different domains 


U, V are both in fi + 


My) 


B 





B 




B 


B 





Di{y) 


B 


B 





Since, for any U(x, •), V(x, ■) E BV L 1 and 7V(f7(-)) + TV(V(-)) + TV{g'{-)) is 
sufficiently small, we have 



C^WUix, •) - < E / h(y)\dy < C \\U(x, •) - V(x, .)|Ui, 

i=l "'- 00 

l<Wi(y)<C , i= 1,2,3,4, 



for some constant Co independent of x and e. Therefore, for any x > 0, 

(4.5) Cf^C/Or, •) - < *(U, V) < C 2 \\U(x, •) - V(x, .)|| L i, 

where Ci and C 2 depend only on 7V(f/(-) + 7V(V'(-) +TF(p'(-)) and the strength 
of the strong shock, which are independent of x. 

Now we examine how the Lyapunov functional $ evolves in the flow direction 
x > 0. Denote Aj the speed of the i— wave <?i(x) (along the Hugoniot curve in the 
phase space) . At a time x which is not the interaction time of the waves either in 
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U or V, 

4 4 

= EE (Wi(v°-Wi(y a -) - \qi(y a +)\Wi(y a +))y a + £ \ qi (b)\W t (b)y b 

a£j i=l i=l 
4 

= {h(y a -)\Wi(y a -)(y a - Xi(y a -)) - W%{y a +)\Wi{y a +){y a - \{y a +))) 

aej i=l 
4 

i=l 

where y Q is the speed of the discontinuity at wave a € J , b = g{x)— stands for the 
points on the boundary, and y b is the slope of the boundary. Define 

(4-6) E a4 = \qt\W+(\+ -y a )-\ q -\Wr(Xr -y a ), 

(4.7) E bii = \qi(b)\Wi(b)(y b - Xi(b)), 

where qf = qi(y a ±), wf = Wi(y a ±), and Xf = Xi(y a ±). Then 

d 4 4 

(4-8) dx~ mix) > V(X)) = £ £ E °" 1 + ^ Eb ' 1 - 

a£j i=l i=l 

Our main goal is to establish the following bounds: 

4 



(4.9) Ei, t i < near the boundary, 

i=l 
4 

(4.10) ^ when a is a strong shock wave in J ', 

i=l 
4 

(4.11) E a , t < 0(l)\a\ when a is a non-physical wave in J ^ 



4 



(4.12) J2 Ea - 1 - °( 1 ) £ l a l when 

a is a weak wave in J . 



i=i 



From (|4.9|l - (|4.12[) . we have 

(4.13) -^$(U(x), V(x)) < 0(l)e. 

ax 

If the constant K2 in the Lyapunov functional is chosen large enough, by the Glimm 
interaction estimates, all weight functions Wi(y) decrease at each time where two 
fronts of U or two fronts of V interact. By the self-similar property of the Riemann 
solutions, $ decreases at this time. Integrating l|4.13|l over interval [0,x], we obtain 



(4.14) ®(U{x), V{x)) < $(E/(0), y(0)) + 0(l)e x. 

In Section 5, we prove (|4.9[) - l|4.12|) . 



WELL-POSEDNESS FOR STEADY SUPERSONIC FLOWS PAST A LIPSCHITZ WEDGE 19 



5. Estimates for the L 1 Stability 

For the case that the weak wave a S J := J{U) U J{V) appears when U and 
V both in £1_ or f2+ and for the case of the non-physical waves in J , estimates 
(|4.11|l - 1|4.12J) can be obtained if \B/(Tq\ is small enough and Ki is large enough, by 
following [Hj. 

Therefore, in this section, we focus the other cases. Cases 1-3 below are all 
related to the strong shock and depend on the wave jump a in U or V; and, by 
carefully adjusting the coefficients Cj and especially relying on estimate 12.15( 1. wc 
can obtain desirable results, which is similar to the Cauchy problem discussed in 
Case 4 is the case near the boundary, which is different from those for the 
Cauchy problem. 

Case 1 (Cross the first strong shock a in U or V): For this case, 

E 1 = BW+(X+ -y a )-\qT\Wr(K -y a ) 

4 

< 0{l)BY,\<k\-KiB\qi\\K-y a \, 



and 

4 4 

^2 Ei = ^(krl(Ar-yc t )(^ + -W i -) + W i +(| g +|(A+- 2 / tt )-| g r|( A r_y Q ))) 

4 4 

< J2 KlB \q+\\X+ -y a \~^52 Kl B\ q r\\\- -y a \. 

i=2 i=2 

Therefore, when n\ is large enough, we have 

4 4 4 

2 



52 Ei <52 Kl B\q+\\X+ -y a \-52 ^ii%r||A- - y a \. 

i=l j=2 i=l 

In Lemma 12.91 let a; = and /3j = 0, which implies <$j = We know that 

4 

Pfc *0(l)Ebrl. * = 2,3,4. 



For the weighted L 1 strength qi(y) in 14.1fl . when q, 1 < i < 4, are larger enough 
relatively to c™, i = 2, 3, 4, we can obtain (|4.1(J|) . 

Case 2 (Cross the weak wave a in between the two strong shocks in U and V): 
For this case, 

£?! = B((W+ -Wr)(\f -y a ) + W?(Xf -Xf)) 
< fl(-«i|a||Ai--y a | + 0(l)|a|). 
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For i = 2,3, 4, 

Ei = \Qt\(Wt-Wn(Xt-ya) + W^(\ q +\(Xf-y a )-\qr\(Xr-y a )) 

< «i|gf||a||Af - y a \ + Kl B((\ q +\ |gr|)(A+ - y a ) + | 9 r|(A+ - K)) 

< Kl | g ±|| £Q ||A±-y Q | +Kl B((|g+|-|gr|)(At-y Q ) + 0(l)|gr||a|). 

Then we have 

4 

$>i < «iO(l)( - \a\ + \a\ J2(kt\ + \%\) + £(1^1 - ICD) + 0(1)|«|. 

Since - II < \q£ - q^\ < 0(l)\a\ when k ^ 1, we can obtain £) i=1 J3j < 
if all the weights cf 1 are sufficiently small and k\ is large enough. 

Case 3 (Cross the second strong shock a in U or V): For this case, by Lemma 
we have 



(5.1) Pi =pi + K siP f. 

Since we have (|2.15l) . the following lemma can be easily obtained. 
Lemma 5.1. There exist c\, c\, and 7 such that 

(5.2) I < 1, 

With Lemma f5. II then we estimate E^: 

Ei = -BWi(\l -y a ) + \q+\W+(\+-y a ) 
< 0(l)B\q+\-KiB\q+\\X+-y a \ 
= 0(l)B\ q +\- Kl Bc a 1 \p+\\X+ -y Q \, 

and, for i = 2,3, 

Ei = -y«){W+ -Wn + W+(\q+\{X+ -y a )-\ q r\{X- -y a )) 

< -niB\q-\(X- -y a ) + 0(l)\q+\ 

< -KiB\qr\(X- -y a ) +0(l)(\qr\ + \q+\). 

By JFHJ and E3, 

Ei = \qi\(X^ -y a )(W+ ~W^) + W+(\q+\(Xi -i Ja )-\q^\(X^ -y a )) 

< KiB (c>4 |(A+ -y a )+ <%K^\p+\(\+ - y a ) - cf \pl\{Xl - y a )) 

< kiB (ct\pi |(A+ - y a ) + 7 c a M\ \*t-V«\- <%\pI\(K ~ V*)) ■ 
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From above, if we choose C4 is small enough relatively to c™ and choose ki is large 
enough, then we obtain 
4 

= -(l-jhiB\qt\\\i-y a \ + 0(l)\q+\ 

i=l 

+n 1 B{cl\p- i \{\+ - y a ) - KA4 - y a )) 

3 

+ Yy^\<k\iK - Va) + 0(1) ■ \q- I) < 0. 

8=2 

Case 4 (near the boundary): For the previous cases, all the desire results depend 
on the cancellation between the two sides of a wave in J . However, it is not the 
case near the boundary since there is only one side near the boundary. Then we 
exploit the exclusive property of the boundary condition l|1.9fl : the flows of U and 
V are tangent to the boundary, which implies that they must be parallel to each 
other. Then we solve the Riemann problem determined by U(b) and V(b). 

Proposition 5.2. Let U(b) — (u,v,p,p) and V(b) — (u,v,p,p) be both in 0(U+), 
v\ju\ — V2/U2 = yb, and v, v ~ 0. Let piib) and Xi be the strength and speed of the 
ith shock in the Riemann problem determined by U(b) and V(b). Then 

(5.4) \y b - Ai| ~|pi(ft)|, i = 2,3, 

(5.5) \ Pl (b)\ < \ P4 (b)\ + 0(1)| P 2(6)||A 2 - y b \ + \ Pl (b)\0(l)\y b \, 

(5.6) |p 4 (6)| = 0(l)|pi(6)|. 

Proof. We divide the proof into two cases. 

Case 1. pi(b) = = p±(b) which corresponds the case p = p: Starting from Ub, 
go along the curves of the second and third families to reach Vf,. These two families 
are the contact Hugoniot curves, and A2 and A3 are constant along the curves. Since 
A2,3 = v /u, r 2 = (1, v/u, 0, 0) T , and r 3 = (0, 0, 0,p) T , v/u keeps unchanged as the 
initial value v(Ub)/u(Ub), i.e. z/f, in this process. Therefore, A2, 3 = yb, i-e., 

2/6 - A 2:3 = 0. 

Case 2. Pxib) 7^ which corresponds to p\ 7^ $4. Starting from U(b), go along 
the first Hugoniot curve to reach Ux, then possibly along the second curve to reach 
U2, the third curve to reach U3, and the fourth Hugoniot curve to reach V(b). 

We project (u,v,p,p) onto the u — v plane to see the relation among pi(b), 
Pi{b), Pz{b), and Pi(b) more clearly. Denote ri| u the projection of ri onto the u 
axis, r2\(u,v) the projection of r 2 onto the u — v plane; and the others are defined 
similarly. At E/+, we have 

r l\u = ~ r 4:\u, r l\v — r i\v, r l|(p : p) = — r i\(p,p)j r 2 = r 2 I (u,^) : r 3\(u,v) = 0. 

The first observation is p±(b) 7^ 0. Since ri|( U)t) ) = k\(— Ai,l) T , the character- 
istic speed is finite and yb ~ 0, so we always have -j^ > y b near U + , i.e. the 



22 GUI-QIANG CHEN TIAN-HONG LI 

derivative dv/du along the 1st curve is always larger than y b in the u — v plane. 
Therefore, v(Ui)/u(U\) ^ v(U b )/u(U b ). On the other hand, we have v(U\)/u(Ux) = 
v(U 2 )/u(U 2 ) = v(U 3 )/u(U 3 ) a,ndv(V b )/u(V b ) = v(U b )/u(U b ). Therefore, v{Ui)/u{U x ) = 
v(JJ'i)/u(JJ'i) — v(Ua)/u(Ua) ^ v(Vb)/u(Vb). To reach V b , there must be some dis- 
tance along the 4th Hugoniot curve. Thus, P4 ^ 0. 

On the u—v plane, we define the signed length of (Ui — Ub) \(u,v) an d (H — £^3)|(u,i>) 
by l\ and 1 4 as follows: 



f \\{Ux-U b )\ {u>v) \\ ifpi>0, 

1 \ - 1| {Ux-Ub) I („,„)|| ifpi <0; 

, _/ ll(H-^3)|(^)|| ifp 4 >0, 

I -||(V6-^3)|( U> „) 



and 

,._/ II(H-^3)|( U ,„) 

„)|| ifp 4 <0. 



The second observation is 

|A 2 -y 6 | = 0(l)|Zi|=0(l)bi(6)|, 
and, since A 2 = v{U\)/u{U\) = v(U2)/u(U2) = A3, we also have 

|A 3 -y 6 |=0(l)|pi(6)|. 

The third observation is 

-k = P2(b) ■ 0(1)(A 2 - y h ) +1 

where lcos9\ — l\cos62, Q\ is the angle between (l,y b ) and r4|( ut ,), O2 is the angle 
between Y\\( UtV \ and (l,y b ), 6\ = 62 + 2/3 for (3 — arctany b , and 

~ cos^ ^ cos(9\ — 2/3) ^ cos6*icos(2/3) + sin9isin{2(3) 

= h (cos(2/3) + 0(l)«n(2/3)) = l x (l + 0(1)/?) = h (l + 0(l)y 6 ) . 
Therefore, we have 

-U = 0(l)p 2 (b)(X 2 - y b ) + h{l + 0(l)y b ). 

Since ri|(„ iP>(!) ) = -r4|( U)PiP ) and n|„ = r 4 | l) at U+, we have 

il — — 
Pi Pa' 

Then we obtain 

(5.7) - Pi (b) = 0(l) P2 (b)(X2 - y b ) +Pi{b){l + 0{l)y b ). 

Therefore, from l|5.7|l . we obtain 

\ P i(b)\ < \ Pi (b)\ + 0(l)b 2 (fo)||A 2 - y b \ + \ Pl (b)\0(l)\y b \, 

and 

|p 4 (&)| = o(i)|M&)|. 

□ 
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By Proposition 15. 21 
Ek,x = \qi(b)\W 1 (b)(y b -\ 1 ) = c a 1 \ Pl (b)\ Kl B(-\ 1 ) + 0(l)\p 1 (b)\ 

< c?|p4(6)|«iB(-Ai)+0(l)|pi(6)|, 
E bti = |ffi(fc)|Wi(6)(^-A i ) = c?|p i (6)|0(l)(i/ 6 -A i ) = 0(l)|pi(6)|, i = 2,3, 
£ M = k 4 (&)|W4(6)(7> fo - A 4 ) = cS|p 4 (6)|KiBAi + 0(l)|p 4 (6)| 

= cSlp^l/siBAx + OWIpiffi)!. 

From Lemma 15 .11 we can find c\ and c\ such that c" < c% . Then, when k\ is large 
enough, we conclude 

4 

J2E b ,i = (cl-cl)\p 4 (b)\K 1 B(-\ 1 )+0(l)\ Pl (b)\ 
i=l 

< (c? - cS) 0(l)bi(6)| KiB(-A x ) + 0(l)bi(6)| < 0. 
6. Semigroup 

We now prove the existence of the semigroup generated by the wave front track- 
ing method. 

Proposition 6.1. If TV(U( )) + TV(g'(-)) is sufficiently small, then, the map 
([/(•),. x) i— > U £ (x,-) := S £ (U(-)) produced by the wave front tracking method is a 
uniformly Lipschitz continuous semigroup with following properties: 

(i) S £ U = U, Sl^ x JJ = S% 1+X2 U; 

(ii) \\SlU-S° x V\\ L x<C\\U-V\\ L x + Cex. 

Proof. Property (i) is obvious since S £ is produced by the wave front tracking 
method. Then we see property (ii). 

Let {U £ } and {V £ } be the front tracking e-approximate solutions of and 
(|1.8[) - (|1.9|) with initial data functions U(-) and V(-), respectively. By 14.14|l and 
()4.5fl . we obtain that, for any x > 0, 

\\U e (x)-V s (x)\\ L i < C'MU £ (x),V £ (x)) 

< Ci$(C/ £ (0),y £ (0)) + CiO(l)ex 

< CiC 2 \\U -V\\ L i + C x O{l)sx. 

This establishes the Lipschitz continuity of the ^-semigroup. □ 

Definition 6.1. Given 5q > 0, define the domain D as the closure of the set 
consisting of the points U : K i— ► K 4 such that there exists one point y l eKso that 
U -U e I^QR.R 4 ) and TV(U -U)< S , where 



U(y) 



U-, y<y\ 

U+, y l < y < boundary. 



24 GUI-QIANG CHEN TIAN-HONG LI 

Remark 6.1. For a solution U(x,y) to the initial-boundary value problem of 
and (|1.8|) - (|1.9[> . if, for any fixed x > 0, U x (y) = U(x,y) G D, then y l = g(0) = 
when x — 0, fcwi y 1 < when x > since f/iere is a strong shock. 

The semigroup denned by the wave front tracking method is set up in the fol- 
lowing theorem. 

Theorem 6.1. IfTV{U{-)) +TV(g'(-)) is sufficiently small, then S £ defined by 
the wave front tracking method is a Cauchy sequence in the L 1 sense. Let S X (U) = 
lim £ ^ S X (U). There exists a constant L such that S : [0, oo) xDnDisa uniformly 
Lipschitz continuous semigroup with following properties: 

(i) SqII — U , S X1 S X2 U — S X1+X2 U ; 

(ii) \\S X U - S X V\\ L , <L\\U-V\\ L i; 

(hi) Each trajectory x i— > S X U yields an entropy solution to the initial-boundary 

problem ifOl) and [fTHj) : 
(iv) If U € D is piecewise constant, then, for x > sufficiently small, the 

function U(x, ■) — S X U coincides with the solution of (|1.3fl and H1.8fl - I|1.9|l 

obtained by piecing together the standard Riemann solutions and the lateral 

Riemann solutions. 

Corollary 6.1. // TV (U (■)) + TV (g' (■)) is sufficiently small, the entropy solution 
to the initial-boundary problem 1)1. 3J) and H1.8(l - (|1.9|l produced by the wave front 
tracking method is unique. 

To prove Theorem 16.11 we need the following lemma which can be found in 6 . 

Lemma 6.2. Let S : [0, oo) x D h- ► D be a globally Lipschitz semigroup. Let X > 0, 
V € D, and V : [0, X] i— > D be a continuous map whose values are piecewise constant 
in the (x, y)— plane, with jumps occurring along finitely many polygonal lines. Let 
L be the Lipschitz constant of the semigroup. Then 
(6.1) 

\\v{x) - s x v\\ L , < l{hv(o) - v\\ L , + [ X ^ mx + h) - ShV(x)hl dx\. 

Proof of Theorem 16.11 It is quite similar to the proof in [H]. The difference 
is that the front tracking method in 6 is to use the cut-off function in the order 
of \fe, while the front tracking method here is to employ the simplified Riemann 
solver rather than the accurate Riemann solver when the interaction term is less 
than £. 

By Lemma f6. 21 
(6-2) 

\\S^V n -S^V m \W < L{\\V n -V m \\ L i+ rimT ll ^ ( ^ K " ) ~^ T/mllLl dx}. 
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Let e m > e n . Then the £ m — approximate solution and e n ~ approximate solution 
only differ when the interaction term of two weak waves or the strength of the wave 
interacting with the strong wave is in [e„, e m \. Suppose that there are N + 1 such 
interactions. For each weak wave interaction between a and f3, 

I a/3 1 = £m, 

and 

either a > y/e^ or (3 > y/e^. 

Let a be large. Then 

S H^(^^) -5(^^)11^ = f; 0( i )em + (i )em 

i=l 
N 

= ^0{l)^a + 0{l)e m 

= 0{l)y/i^TV(V(-))=0(l)y/^. 

Therefore, S x n V n is a Cauchy sequence, which converges in the L 1 sense. Hence, 
the map 5 : [0, 00) x D h B as the limit of the approximate solutions produced by 
the wave front tracking method is well-defined. 

Next, we prove (i) to (iv). Facts (i), (ii), and (iv) are obvious since S is the 
limit of S e produced by the wave front tracking method. It is similar to prove 
(iii) as the only difference is that the wave front tracking method we employ 
here is slightly different. Finally, we can see that the entropy solution satisfies the 
boundary condition due to the construction of our approximate solutions. This 
completes the proof of Theorem 6.1. 

7. Uniqueness of Entropy solutions in a broader class 

In this section, we first prove the semigroup S defined by the wave front tracking 
method is the only standard Riemann semigroup (SRS) which is defined as Def- 
inition |7TT] That is, the semigroup defined by the wave front tracking method is 
the canonical trajectory of the standard Riemann semigroup (SRS). Then we prove 
that the uniqueness of entropy solutions in a broader class, i.e. the class of viscosity 
solutions as defined in 0]. The main point is to prove that, in the class of viscos- 
ity solutions, the entropy solution is unique, which coincides with the trajectory 
produced by the wave front tracking method. 

Definition 7.1. We say that problem (jl.3JI and 1)1.8(1 - (|1.9[1 admits a standard 
Riemann semigroup (SRS) if, for some 5q, there exists a continuous mapping: 
R : [0, 00) xBhB and a constant L with the following properties: 

(i) RqU = U, R Xl R X2 U — R Xl+X2 U : 

(ii) \\R X U - R x V\\ L i <L\\U-V\\ L i; 
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(iii) If U G H> is piecewise constant, then, for x > sufficiently small, the 
function U{x, ■) = R X U coincides with the solution of 1(1.3(1 and 1(1.8(1 ~ l(1.9() 
obtained by piecing together the standard Riemann solutions and the lateral 
Riemann solutions. 

Theorem 7.1. Let problem 1(1.3(1 and ((1.8(1 - 1(1.9(1 admits a standard Riemann semi- 
group R : [0, oo) x B i— > D. Let S be the semigroup generated by the wave front 
tracking method, i.e. S X {U) = lim e ^o S X {U). IfU^zO, then R X U — S X U for all 
x > 0. 

The proof of the theorem is similar to the proof in [3] by using Lemma l6.2l and the 
fact that, locally in x direction, the wave front tracking method and the standard 
Riemann semigroup (SRS) both have the structure of the Riemann solutions. 

As in 0], there are two types of local approximate parametrices for (|1.3ll : One 
is derived from the self-similar solution of a Riemann problem, and the other is 
obtained by "freezing" the coefficients of the corresponding quasilinear hyperbolic 
system in a neighborhood of a given point. 

Let U : [0, oo) x R i— > R 4 be a function. Fix any point (r, £) in the domain of U. 
If U(t, •) € B, then the bound on the total variation implies the existence of the 
limits 

U~ = lim U(T,y), U + = lim U(r,y). 

Denote by to = u>(x,y) the corresponding solution of the Riemann problem with 
U~ and U + and by A a upper bound for all characteristic speeds, i.e., 

(7.1) sup|Ai([/)| < A, i = 1,2,3,4. 

u 

For x > t, define the function 

(7 9) W* It ,A - / w(aj-r,y-0 if|y-f| < A(x-t), 
( J W wV>V)-[ U(T>y) i{\y-Z\>\( X -T). 

Set A=DW(U(t, ^)) and B^DH(U(t,^)) the Jacobian matrices computed at 
the point U(t, £). For x > r, define T / F^ (7t ^- ) as the solution of the linear Cauchy 
problem with constant coefficients 

(7.3) AV x + BV y =0, V(T,y) = U(r,y). 

Then the functions W# and W h depend on the values U(t,£) and U(t, £±). Next, 
we define viscosity solutions which have the same local characterization as W# and 
W b . 



Definition 7.2. Let U : [0,AT] i— ► B be continuous with respect to the L 1 norm. 
We say that U is a viscosity solution of system p. 3(1 and 1)1.8(1 - 1(1.9(1 if there exist 
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constants C and A satisfying (17.1(1 such that, at each point (r, Z) £ [0, X) x M, iw/ien 
p and e are sufficiently small, 



\U(t + e,y) -W* ri) {x,y)\dx 



(7.4) < CTV{U(r) : (Z-p,Z)^(Z,Z + p)}, 



1 



\U{r + e,y)-W b {UT0 {x,y)\dx 

(7.5) < C(TV{U(r) : (£ - p, Z + p)}f . 

Lemma 7.1. Lef (a, 6) 6e a (possibly unbounded) open interval and let R be the 
standard Riemann semigroup (SRS). If U , V € D, then, for all x > 0, 

(7.6) I \R x U(y)-R x V(y)\dy<L f \U (y) - V (y)\dy . 

Proof. If two initial data functions U and U € B coincide on (a, 6), due to the finite 
speed of propagation, then the semigroup generated by the wave front tracking 
method has S x U(y) = S x U(y) for y £ (a + Xx,b~Xx). By Theorem 177X1 R x U{y) = 
S x ll{y) = S x U{y) = R x U{y) for all y e {a + \x,b-Xx). Next, for U,V £B, define 

U = U X[a,b] + U-X(-oo,a) + U+X(b,g(x)) ) 

where x is t ne characteristic function, and U- and [/+ are the states in the definition 
of B. It follows that U G B. Similarly, V = VX[ a ,b) + U-X(-oo,a) + U + X(b, g {x)) 6 B. 
The uniform Lipschitz continuity of the semigroup R implies 

b—\ X r-b—\x 

\R x U(y)-R x V(y)\dy = / \R x U(y) - R x V(y)\dy 

a-\-Xx Ja-\-Xx 

< \\R x U{y)-R x V{y)\\ 

b 



< L\\U-V\\=L f \U~V\dy 

J a 



□ 



Theorem 7.2. Assume that problem (|1.3|) and l|1.8fl - l|1.9|) admits a standard Rie- 
mann semigroup R. Then a continuous map U : [0, X] i— > B is a viscosity solution 
of i|T3|) and lfT"%j) - ifOl) z/ and on/y «/ 

(7.7) U{x,-) = R X U for any x e[0,X}. 

Corollary 7.1. For system l|1.3fl and l|1.8fl ? the entropy solution is unique in 
the class of the viscosity solutions, which coincides with the trajectory S X U generated 
by the wave front tracking method, i.e. a continuous map U : [0, X] i— > B is a 
viscosity solution if and only if 

(7.8) U(x, •) = S X U for any x e [0, X}. 
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The proof is similar to the argument in 0]. The only difference is that there is 
a strong shock in our case; however, we can still carry out the proof as long as the 
convergence of the wave front tracking method is achieved which has been proved 
in Section 3. 

Remark 7.1. For the potential flow, isentropic or isothermal Euler flow (|1.6|l . 
which are the simpler cases as the L 1 stability problem as concerned, we obtain the 
same results as the full Euler equations. 
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